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Binary systems of massive black holes will be detectable by the Laser Interferometer Space An- 
tenna (LISA) throughout the entire Universe. Observations of gravitational waves from this class of 
sources will have important repercussions on our understanding of the behaviour of gravity in the 
highly non-linear relativistic regime, the distribution and interaction of massive black holes at high 
redshift and the formation and evolution of cosmic structures. It is therefore important to address 
how accurately LISA can measure the source parameters and explore the implications for astronomy 
and cosmology. 

Present observations and theoretical models suggest that massive black holes could be spinning, 
possibly rapidly in some cases. In binary systems, the relativistic spin-orbit interaction causes the 
orbital plane to precess in space producing a characteristic signature on the emitted gravitational 
waves. In this paper we investigate the effect of spins on the gravitational wave signal registered at 
the LISA output - we provide ready-to-use analytical expressions of the measured signal - and the 
implications for parameter estimation. We consider the in-spiral phase of binary systems in circular 
orbit undergoing the so-called "simple precession" and we approximate the gravitational radiation 
at the restricted post'^'^-Newtonian order. We show that the presence of spins changes dramatically 
the signature of the signal recorded by LISA. As a consequence, the mean square errors associated 
to the parameter measurements are significantly smaller than the ones obtained when the effect of 
spins is neglected. For a binary system of two 10^ Mq black holes, the angular resolution and the 
relative error on the luminosity distance improve by a factor « 3-to-lO; the fractional errors on the 
chirp mass and the reduced mass decrease by a factor 10 and ~ 10^, respectively. 



I. INTRODUCTION 



Binary systems of massive black holes in the mass range 10^ Mq — 10^ Mq are among the most spectacular sources 
that will be observed by the Laser Interferometer Space Antenna (LISA), an ES A/NASA space-borne laser interfer- 
ometer aimed at observations of gravitational waves (GW's) in the low-frequency region of the spectrum ~ 10~^ Hz - 
0.1 Hz, and expected to be launched in 2011 y|. The interest in the detection of gravitational radiation from massive 
black holes (MBH's) is two-fold 0: (i) for fundamental physics, because it will provide a thorough understanding 
of the behaviour of strong gravity in the highly non-linear relativistic regime, and (ii) for astronomy and cosmology 
because GW's provide an "orthogonal" and complementary mean, with respect to conventional astronomy, of inves- 
tigating the high-redshift formation and evolution of cosmic structures, such as galaxy interactions and mergers, and 
exploring the demography of massive black holes. 

Massive black holes (MBH's) are considered to be essentially ubiquitous in the centre of galaxies in the local and 
low-redshift Universe Our Galaxy and NGC 4258 0] provide probably the most solid observational evidence 

for the existence of such black holes Q ■ The evidence for the presence of black holes in galaxies at high redshift is 
more circumstantial, but the fact that galaxies, at some stage of their life, harbour AGN's, and we observe quasars 
at z ^ 6 strongly favour the hypothesis that black holes arc fairly common also at high redshift. This conclusion 
is supported by recent results ^ that show that the high-redshift quasar distribution and the present day density 
of black holes can be reconciled by invoking galaxy mergers where the relation between host velocity dispersion and 
black hole mass is consistent with the one observed locally ^3- The present day population of MBH's is therefore 
likely to be the dormant remnant of those earlier activities. Two black holes are brought together at the centre of a 
merger remnant by dynamical friction from the stars of the common host. The further evolution of a binary must 
then take place on a time scale short enough for the system to coalesce within an Hubble time due to radiation 
reaction. Two are the main mechanisms that could make the binary hard enough to satisfy this condition: dynamical 
friction from stars in the core 0, 0, 0, [l4|, [l5| and gas dynamical effects 16, 17]. However, despite all these 
circumstantial evidences, the rate of coalescence of massive black hole binaries in the Universe is still controversial: 
present astrophysical estimations vary widely from ~ 0.1 yr~^ to ~ 100 yr^-^ 0, E, 113, |^ [12, [13, 01 1 depending 
on model assumptions and BH mass range. At present there are no direct observational evidences for the presence of 
massive black holes binary systems, with the possible exception of OJ287 [25| . 

LISA will be able to detect MBH binaries at a moderate-to-high signal-to-noise ratio throughout the entire Universe 
(redshift z « 10 and beyond, if black holes are already present). During the whole LISA observational campaign - 
the nominal duration of the mission is 3 years, but the instrument is expected to operate for about 10 years, if no 
catastrophic failure occurs on the spacecraft - one can therefore expect from a handful to a very large number of 
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detections. 

LISA is an all-sky monitor. At any one time, the instrument maps the whole sky, although with a different 
sensitivity depending on the location of the source and the polarisation of the waves. By coherently tracking the 
phase and amplitude evolution of GW's, one can extract precise information about a source. In observations of binary 
systems, the number of parameters on which the waveform depends, and that one needs to extract from the data, is 
large, in the most general case 17. Clearly correlations among parameters are inevitable, degrading the precision of 
the measurements. 

A few studies have been carried out so far to address how accurately LISA can measure the source parameters and 
the implications for astronomy and cosmology. Because of the complexity of the problem and open theoretical issues 
in the modelling of the gravitational waveform in the most general case, assumptions, possibly ad-hoc, have been 
introduced on the waveform model. Cutler |26l | has considered the in-spiral phase of binary systems approximating 
the waveform at the restricted post-'^-^-Newtonian order. Sintes and collaborators ll^,!!^ ^^'^ Moore and Hellings |29| 
have investigated the parameter estimation including post^-Newtonian corrections to the amplitude and the phase. 
Hughes [sSl has considered restricted post^-Newtonian waveforms for the in-spiral phase and has included the ring- 
down portion of the signal; more recently Seto has included the effect of the finite length of the LISA arms using 
restricted post^'^-Newtonian waveforms |31j . All these analysis have a priori assumed that black holes have no spins, or 
the spins are parallel to each other and the orbital angular momentum. However, present astrophysical observations 
and theoretical models suggest that massive black holes are likely to be rotating, possibly rapidly in some cases. 
Observations of AGN's and quasars are highly consistent with rapidly rotating black holes; observations of black holes 
in quiescence are in agreement with a slow, if any, rotation and some theoretical models suggest small-to-moderate 
spins (see, e.g., [s^l)- If there is little doubt that black holes are, in general, rotating, how strong spins are is to a 
large extent unknown. Quite likely, gravitational wave observations shall provide the most effective tool in addressing 
this issue. Spins are likely to add vast complexity to the signal recorded at the detector output, and require more 
parameters to describe the waveform [s^, 13 HE HEl ■ In this paper we show that spins do change dramatically the 
structure of the signal observed by LISA and, as a consequence, affect significantly (even by orders of magnitudes for 
some selected parameters) the errors with which LISA can measure the source parameters. 

In this paper we restrict our attention to the in-spiral phase of massive binary systems of comparable mass. We 
consider circular orbits, model the waveform at the restricted post^'^-Newtonian order and assume that binary systems 
undergo the so-called simple precession [s^. which is the relevant scenario from an astrophysical viewpoint. The goal 
of the paper is two fold: (i) to provide ready to use analytical expressions for the LISA detector output when binary 
systems undergo spin-orbit precession, and (ii) to show the dramatic change in parameter estimation that occurs when 
spin-orbit modulation in phase and amplitude is taken into account. The actual parameter space is so vast - the 
waveform is described by 12 parameters - that here we concentrate on typical systems of two 10^ Mq black holes. A 
thorough exploration of the whole parameter space requires large-scale, CPU-intensive Monte-Carlo simulations and 
is beyond the scope of this paper; such analysis is currently in progress and will be reported in a separate paper |37| . 

The organisation of the paper is as follows. In section^ we review the main features of the LISA mission, in 
particular its orbital configuration which is central for the topic discussed in this paper. In Section [llll we review the 
main properties and key equations regarding the emission of GW's from in-spiralling binary systems, with emphasis 
on the role plaid by spins, and we stress the differences with respect to the non-spinning, i.e. non-precessing, case. 
Section IIVI and contain the key results of the paper. In Section IIVI we derive explicit ready-to-use analytical 
expressions for the signal measured at the LISA detector output for binary systems in circular orbit that undergo 
simple precession; the waveform is computed within the restricted post^'^-Newtonian approximation. In section Ivl we 
explore the errors with which the twelve parameters of the signal can be determined and compare the results with 
the case where spins are neglected. In Section IVTl we summarise our conclusions and present pointers to future work. 

II. THE LISA MISSION 

In this section we review the main properties of the LISA mission. We refer the reader to for more details. 

A. The LISA orbit 

LISA is an all-sky monitor with a quadrupolar antenna pattern. Its orbital configuration is conceived in order to 
keep the geometry of the interferometer as stable as possible during the mission, which in turn provides a thorough 
coverage of the whole sky: a constellation of three drag- free spacecraft (containing the " free- falling test masses") is 
placed at the vertexes of an ideal equilateral triangle with sides ~ 5 x 10^ km; it forms a three-arms interferometer, 
with a 60° angle between two adjacent laser beams. The LISA orbital motion is as follows: the barycentre of the 
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instrument follows an essentially circular heliocentric orbit, 20° behind the Earth; the detector plane is tilted by 60° 
with respect to the Ecliptic and the instrument counter-rotates around the normal to the detector plane with the 
same period lyr. 

Following ^26\ wc introduce two Cartesian reference frames "s^: (i) a fixed " barycentric" frame {x,y, z) tied to the 
Ecliptic and centred in the Solar System Barycentre, with z perpendicular to the Ecliptic, and the plane {x, y) in 
the Ecliptic itself; (ii) a detector reference frame {x',y',z'), centred in the LISA centre of mass and attached to the 
detector, with z' perpendicular to the plane defined by the three arms and the x' and y' axis defined so that the unit 
vectors Ij {j = 1,2,3) along each arm read 



Ij — cos 



TT TT , 

T2 + 3(^'-^^ 



X + sm 



T2 + 3(^-^) 



(1) 



In the Ecliptic frame the motion of LISA's centre-of-mass is described by the polar angles 



where 
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lyr 



(2) 



(3) 



and $0 sets the position of the detector at same arbitrary reference time. The normal to the detector plane z' precesses 
around z according to 
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cos $(i)x + sin cE>(i)y] 
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The time evolution of the unit vectors Ij {j = 1,2,3) along each arm is described by the following expression [2 



(4) 



1, 



- sin a j{t) cos^{t) — cos a j{t) sin$(t) 



- smaj{t) sin$(i) + cos aj{t) cos3>(t) 



sma_,(r) 



where aj (i) increases linearly with time, according to 

a^(t) = riei - (j - l)7r/3 + olq , 
and ao is just a constant specifying the orientation of the arms at the arbitrary reference time i = 0. 



(5) 



(6) 



B. The LISA detector output 

The strain h(€) produced at the output of the LISA Michelson interferometer by a GW signal characterised by the 
two independent polarisation states /i+(i) and h-^if) is 

^ ^ [F|^)(t)/,+ (t) +4')(i)/ix W] . (7) 

In Eq. lO F^_^ (t) and F^^\t) are the time dependent antenna patters and the factor \/3/2 ~ sin(7r/3) comes from the 
60° opening angle of the LISA arms. and vary with time because during the observation the interferometer 
changes orientation with respect to the source; in fact, in-spiral binaries are long-lived sources in the low frequency 
band. The index l ^ 1 ,11 labels the two independent Michelson outputs that can be constructed from the readouts 
of the three arms if the noise is uncorrelated and " totally symmetric" p3 | . They are equivalent to the outputs of two 
identical interferometers in the same location, rotated by 45° one with respect to the other. 

The functions describing the interferometer beam pattern depend on the source location in the sky N (GW's 
propagate in the — N direction) and wave polarisation, which is related to the orientation of the orbital plane, and 
therefore the orbital angular momentum L, with respect to the detector. With respect to the frame tied to the Solar 
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system barycentre, N and L are described by the polar angles {9n jiPn) and {0^,(1)^). Equivalently, in the frame 
attached to the detector, N and L are identified by the angles , (p'j^) and (0^ ^^'^ reference frame attached 

to LISA, the antenna beam patters read 

F+ {e'j^ , (/)'^ Vjv) = ^(1 + cos 6*^) cos 20^ cos2V'^ - cos 6*^ sin 20'^ sin2V'^ , (8a) 

Px [O'nA'n^'n) = ^(l + cos6'^)cos2<?!)^sin2V'5v + cos6'5vsin2<?!)'^cos2?A5v, (8b) 

where V'jv describes the "polarisation angle" of the waveform in the detector frame. The antenna patterns H8a() 
and H8b|l for the two outputs t = J and l = II are therefore 



F+^^{9'j,,{4>'j,^7t/A) {i = II) 



F^Ut) = { . (9) 



Note that because of the detector change of orientation, the angles 0^, 0'^ and V'jv are time dependent. As a function 
of the angles measured in the Solar System Barycentre they read 

1 

cos9'j^{t) = -cos^Tv — sinON cos{^{t) — , (10) 

„ TT ,1 [ ^cos0jv + ^singjvsin(^(t)-,/>jv) | 

12 I smc/jv sm($(i) — 0jv) I 

,, L • z' - (L • N) (z' • N) 
N • (L X z') 



where 



= Tiei- ^ - ^(j - 1) + So, (13) 



and So sets the orientation of Ij at i = 0. We are not yet ready to derive the explicit expression for ip'j^, Ea. (|12|l . 
because for spinning black holes, the source at the centre of this paper, L is not a constant of motion. We derive the 
necessary equations in the following section. 

III. GRAVITATIONAL WAVES FROM BINARY SYSTEMS 

In this section we briefly review, partly to fix notation, the basic concepts and formulae regarding the emission of 
gravitational waves by binary systems, with emphasis on the effects produced by spins. We refer the reader to [ssLIb^ 
for details and to 39] and references therein for a thorough review regarding post-Newtonian gravitational waveforms. 

The whole coalescence of a binary system is usually divided in three distinct phases 40]: (i) the adiabatic in-spiral, 
during which the BH orbital evolution occurs on the time-scale of the gravitational radiation reaction, which is much 
longer than the orbital period; (ii) the merger, which takes place at the end of the in-spiral when the binary orbit 
becomes relativistically unstable and the black holes enter a free fall plunge; and (iii) the ring-down, when the resulting 
BH settles down in its final stationary Kerr state, oscillating according to its normal modes. In this paper we consider 
only the in-spiral phase. In fact the time to coalescence for a system radiating at frequency / is: 



1.2x10^ / ^-/VM(l + .)^-/^ 



10-4 Hz/ V 10^ Mg 



where M and rj are the total mass and the symmetric mass ratio, respectively, defined in Eqs. (|16b|l - H16e|l . r is the 
observed life-time of a binary system as recorded in the solar system for a source at redshift z with intrinsic total mass 
M that enters the detector window at the observed frequency / = 10^* Hz. Binary systems are therefore long-lived in 
the LISA band, and one critically relies on long integration times to disentan gle t he sources parameters, in particular 
to resolve the source position in the sky and measure its luminosity distance |26| . 
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As we consider only the in-spiral phase of the coalescence, in the frequency domain we shut-off the signal at 
frequency 

f. 1 

J ISCO 



7r63/2M(l-f z) 



3 fM{l + z) 



-1 



which corresponds to the innermost stable circular orbit of a particle orbiting a Schwarzichild black hole (ISCO). The 
real transition from in-spiral to merger will occur at a frequency somewhat (by a factor ^ 2) different from /isco, but 
due to the difficulty of defining the ISCO we will adopt the value (|15|l . The results presented in this paper are anyway 
essentially unaffected by changing the cut-off frequency by a small factor. 

We also assume the orbit to be circular. This is completely reasonable from an astrophysical point of view, as 
a massive binary system of roughly equal mass black holes has enough time - regardless of the initial eccentricity 
at which it forms - to circularise before it enters the observational window, due to dynamical friction and radiation 
reaction [Tllilll. 



A. The amplitude and phase evolution 



We consider two massive black holes of mass mi and m2, and spins Si and S2, respectively. Several mass parameters 
are actually useful: 



M 
M 



mi+m2, 

mi 777,2 

M ' 
^3/5^,^2/5 ^ ^3/5 

= ■ 

M ' 



(16a) 
(16b) 

(16c) 
(16d) 
(16e) 



they represent, respectively, the total mass, the reduced mass, the chirp mass and the symmetric mass ratio. 

The system evolves by loosing energy and angular momentum through emission of gravitational waves of increasing 
frequency and amplitude. The emitted radiation is given by the super pos ition of harmonics at multiples of the orbital 



period, and the two independent polarisation amplitudes and |42l | can be schematically represented as |4 



(17) 



where are time dependent functions of the source parameters and (/)oib(0 the binary orbital phase. The 

strongest harmonic is associated to the quadrupole moment of the source, and therefore corresponds to fc = 2. In 
this paper we consider the so-called restricted post-Newtonian approximation to the radiation (|17|l : we retain only 
the leading order contribution to the amplitude, therefore only H]_ ^ at the Newtonian order, and take into account 

'orb(i)- In this approximation Eq. (|17|l reads: 

(18a) 



post-Newtonian corrections only to the phase of the signal (f>(t) — 2 

A^5/3 



hx = -4 



Dl 

^5/3 



1 + L • N 



[nff^ cos<^(i) 



Dr 



L-N) (nff sin <f>{t), 



(18b) 



where is the source luminosity distance j4J| , N the unit vector pointing toward the binary centre of mass (in our 
convention the GW propagation direction is — N) and L the orbital angular momentum unit vector. 

We review now the expression for the GW phase (/'(t). The signal frequency /, twice the orbital one, evolves 
according to 



!jL = ^7rS/3_A/(5/3 Hl/3 

dt 5 

/34103 13661 
^ V 18144 ^ 2016 



1 - + ) i^Mff^' + (47r - (3)inMf) 

59 



(19) 
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where Eq. (|19() is valid through the post^-Newtonian order and 

i=l L 



L 



48 



-247 



Si 



721 



L 



Si 



L 



Si 



(20) 
(21) 



are the so called spin-orbit and spin-spin parameters, respectively 45]. From Eq. (|19|l and d(j)/dt = 2i:f, one can 
derive the time and phase evolution of the gravitational radiation: 



t{f) = ic-5(87r/)-8/3^-5/3 



3 V336 4 
/3058673 5429 617 



n4/3 



V 1016064 1008 ' 144 
_ /3058673 5429 617 



(22) 



V 1016064 1008 



N4/3 



(23) 



tc and 4>c, the time and phase at coalescence, are constants of integration, defined as the value of t and (j) at (formally) 
f — oo. As we mentioned, massive black holes spend months-to-years in the LISA observational window, depending 
on the mass. For future reference, if we assume that LISA observes the final year of in-spiral of a binary, then the 
signal sweeps the frequency band between 



/a ^ 4.2 X 10 



-5 



To 



bs 



lyr 



-3/8 



M{l + z) 



106 Mq 



-5/8 



(24) 



which we call "arrival" (or "initial") frequency, and /isco- An interesting quantity to consider is the total number of 
gravitational wave cycles recorded by the detector: 



Afit)^^ cj,{t')dt'. 

From Eq. H23|l and (|24|l . it is easy to compute the number of cycles in the final year of in-spiral: 

\-\ -5/3 



(25) 



w 4.8 X 10^ 



2.1 X 10-^ 



fa 



10-4 Hz 



^-5/3 


'Mil + z)' 




106 Mq 



^5/8 


■^1(1 + ^)" 




106 Mq 



-5/8 



(26) 



Table ^ contains reference values for /a, /isco and Af for a number of choices of masses and observation times. 



B. Spin induced precession 



Present observations and theoretical prejudices suggest that massive black holes could be rapidly rotating. As we 
have mentioned in the previous section, if black holes are spinning then the GW phase H23() contains two additional 
terms, one at the 1.5PN order, which is proportional to p, and the other at the 2PN order, which is proportional 
to a. However spins introduce a much more dramatic effect in the structure of the waveform, which is qualitatively 
different from the ones just mentioned: the orbital angular momentum L and the spins Si and S2 change orientation 
due to spin-orbit interaction (1.5PN effect), so that the radiation detected has a time varying polarisation. L, Si and 
S2 precess, over a time scale longer than the orbital period (2// ^ 30 min) but shorter than the observation time 
(^ 1 yr), around the (almost) fixed direction of the total angular momentum L-I-Si -I-S2. This effect is actually quite 
dramatic and drastically changes the signature of the detected signal, as the orientation of the orbital plane changes 
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TABLE I: Summary of the frequency band swept by massive black hole binary systems and the number of wave cycles recorded 
at the detector output in LISA observations for selected masses and integration times. The table shows the initial and final 
frequency of the GW's, /a and /isco respectively, the corresponding time of observation Tobs, and the number of waves cycles 
recored at the output of the detector, coming from the various contributions of the terms in Eq. 1231 : Newtonian: the first 
term; post ^-Newtonian (IPN): the second term, proportional to {tt M f)^^^; tail: — 10 vr (tt A/ /); spin-orbit: (5/2) /3 (tt M /); 
post^-Newtonian (2PN): the last term, proportional to (ttM with a = 0; spin-spin: the -5cr(7rAf Z)*/^ term. The 
sources are selected with masses in the ranee 10'^ Mq - 10^ Mq- for each mass choice, three different /a's are considered: the 
first corresponds to Tobs =1 yr, which does depend on mi and m2, whereas the other two are the same for all mi and m2, and 
correspond to /a = 5 x 10~^ Hz , and 10~* Hz, respectively. Notice that depending on /a, mi and m2, systems are observable 
for different times. The fiducial source is at 2: = 1. 















number of wave 


cycles 




mi 


m2 


fa 


/isco 


Tobs 


Newt. 


IPN 


tail 


spin-orbit 2PN 


spin-spin 


(Mq) 




(xlO"*Hz) 


(xlO"*Hz) 


(yr) 








(x/3) 




(xa) 


10^ 


10^ 


0.070 


1.099 


1.000 


342 


42 


-40 


3 


3 


-1 


10^ 


10^ 


0.500 


1.099 


0.004 


10 


3 


-5 





1 





10^ 


10® 


0.153 


1.999 


1.000 


761 


87 


-103 


8 


7 


-2 


10^ 


10" 


0.500 


1.999 


0.038 


96 


22 


-34 


3 


3 


-1 


10^ 


10" 


1.000 


1.999 


0.005 


23 


7 


-13 


1 


1 





10^ 


10^ 


0.350 


2.177 


1.000 


1792 


290 


-455 


36 


32 


-10 


10^ 


10^ 


0.500 


2.177 


0.366 


949 


187 


-319 


25 


24 


-8 


10^ 


10^ 


1.000 


2.177 


0.045 


238 


65 


-130 


10 


11 


-4 


10^ 


lO"* 


0.776 


2.196 


1.000 


4102 


990 


-1881 
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153 


-51 


10^ 


10* 


0.500 


2.196 


3.649 


9477 


1834 


-3160 


252 


235 


-78 


10^ 


10* 


1.000 


2.196 


0.447 


2383 


647 


-1296 


103 


110 


-37 


10'' 


10" 


0.295 


10.993 


1.000 


1453 


102 


-76 


6 


5 


-1 


lO'^ 


10" 


0.500 


10.993 


0.246 


602 


59 


-51 


4 


4 


-1 


lO'^ 


10" 


1.000 


10.993 


0.038 


187 


28 


-30 


2 


3 


-1 


lO'' 


10^ 


0.647 


19.987 


1.000 


3216 


215 


-201 


16 


11 


-3 


10'' 


10^ 


0.500 


19.987 


1.991 


4950 


280 


-243 


19 


12 


-4 


10" 


10^ 


1.000 


19.987 


0.312 


1552 


137 


-144 


11 


9 


-3 


10" 


10* 


1.508 


21.768 


1.000 


7571 


748 


-943 


75 


55 


-18 


10® 


10* 


0.500 


21.768 


19.285 


48121 


2366 


-2176 


173 


96 


-31 


10" 


10* 


1.000 


21.768 


3.019 


15096 


1155 


-1300 


103 


68 


-22 


10^ 


10^ 


1.242 


109.929 


1.000 


6161 


247 


-141 


11 


8 


-2 


10^ 


10^ 


0.500 


109.929 


11.261 


28097 


618 


-266 


21 


11 


-2 


10^ 


10^ 


1.000 


109.929 


1.781 


8848 


308 


-165 


13 


9 


-2 


10^ 


10* 


2.729 


199.872 


1.000 


13594 


519 


-374 


30 


16 


-5 


10^ 


10* 


0.500 


199.872 


91.802 


230221 


2867 


-1208 


96 


33 


-9 


10^ 


10* 


1.000 


199.872 


14.494 


72508 


1430 


-753 


60 


25 
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during the year long observation time, and the angular momenta complete several precession cycles around the fixed 
detection. 

The spin-orbit induced precession is the key new physical effect which is introduced in this paper. Here we review 
the main concepts and expressions and refer the reader to |3,li l34| for more details. In this paper we restrict our 
attention to the so-called simple precession ,33j , which takes place when (i) the BH masses are equal (mi W2), 
or (ii) one of the BH's has negligible spin (for convention, in this paper we shall assume 5*2 = 0). Furthermore, 
one needs to impose the condition that the angular momenta are oriented so that L 7^ —(Si + S2). Under this 
assumptions the equations that describe the evolution of the angular momenta simplify considerably, and one can 
construct explicit analytical solutions for the relevant quantities. Moreover the regime that we consider is well justified 
form an astrophysical point of view. 

For circular orbits through the post^'^-Newtonian order the equations describing the evolution of L, Si and S2 read 



8 



3 ^2 + 1^) J 
3 \ 2mi / 



X L . 



5* = 



3m2\ 



— ( 2 + I 



X S, 



dt 



(Si -82) =0, 



(27) 



(28) 



(29) 



S = 0, (30) 

where 

J = L + S (8 = 81 + 82) (31) 



is the conserved total angular momentum |46l l47l | (indeed L = —8). 

Following [33 |. we highlight the main features of simple precession and provide, in a ready-to-use form, the relevant 
quantities needed for the computation of the LISA detector output. We will neglect the post^-Newtonian corrections, 
while retaining the post-'^'^-Newtonian terms (spin-orbit), which are dominant, for the mi — m2 case; they vanish 
anyway for 5*2 = IH, |3 ■ 

One can summarise the main features of the simple precession as follows. J, L and 8 precess with the same angular 
velocity 

".-f^ + S^U (32) 



around the fixed direction 



where 



The quantities 81-82, 5* = [81 + 82 1 and 



2toi 



Jo = J - eJ X L , (33) 



8 • L (35) 



are constant during the in-spiral. The phase of the gravitational signal contains an additional contribution, the 
so-called " Thomas precession phase" , given by (13 

5p(j){t)^- j^' 5p^{t')dt' , (36) 

where 

5Mt) = ^l^. {LxN)-L. (37) 

In analogy with the number of wave cycles, we can introduce here the number of precession cycles Afp that L and 
8 undergo around J from some initial frequency fa, at time t^, to the final coalescence at tc- 

1 



Afp = — npit)dt. (38) 

Z TT ' 
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TABLE II; The total number of precession cycles Mp and the precession angular frequency 57p at the beginning and end of 
LISA observations, for the final year of in-spiral of a binary system at redshift z = 1 with k = S ■ L = 0.9 for selected values of 
the masses and spin magnitude. 



7Tl\ 
(Mr,') 


7722 




at L 
(xlO"" Hz) 


0„ at 

J t-p CXiVj J ISCO 

(xlO'^Hz) 


Number of precession cycles 


10'^ 


10' 


0.95 


0.009 


1.254 


11 


10'^ 


10' 


0.50 


0.007 


0.892 


g 


10'^ 


10'' 


0.10 


0.005 


0.579 


7 


10'^ 


10*' 


0.95 


0.008 


1.043 


11 


10'^ 


10® 


0.50 


0.005 


0.632 


7 


10'^ 


lO'' 


0.10 


0.003 


0.272 


4 




1 






n qi 1 




10'^ 


10^ 


0.50 


0.014 


0.489 


19 


10'^ 


10^ 


0.10 


0.003 


0.114 


5 


10'^ 


10*1 


0.95 


0.133 


0.895 


127 


10'^ 


10** 


0.50 


0.064 


0.472 


71 


10'^ 


10** 


0.10 


0.012 


0.096 


15 


lO" 


lO*" 


0.95 


0.018 


12.540 


25 


lO" 


lO'' 


0.50 


0.015 


8.919 


20 


lO" 


lO'' 


0.10 


0.013 


5.788 


16 


lO" 


10^ 


0.95 


0.015 


10.432 


23 


10® 


10^ 


0.50 


0.011 


6.319 


16 


lO" 


10^ 


0.10 


0.007 


2.716 


g 


lO" 


10** 


0.95 


0.047 


9.110 


74 


lO" 


10** 


0.50 


0.025 


4.892 


40 


10" 


10^ 


0.10 


0.007 


1.144 


11 


10^ 


10^ 


0.95 


0.038 


125.396 


54 


10^ 


10^ 


0.50 


0.034 


89.194 


46 


10^ 


10" 


0.10 


0.030 


57.877 


39 


10^ 


10^ 


0.95 


0.030 


104.315 


48 


lO'^ 


10^ 


0.50 


0.023 


63.186 


34 


lO'-' 


10^ 


0.10 


0.016 


27.159 


22 



Depending on the mass ratio of a binary system, whether mi ~ r7i2 (L ^ S) or mi ^ m,2 {L <C 5*), A/'p is "well 
approximated by the following expressions: 

■I^P ~ \ , \ / \ / r. \ / \ -2/3 / , \ -2/3 • l^JJ; 

i 2 (l + If^) (^) {i.) (to^) (to^) [L « S) 

Notice that the precession cycles accumulate at low frequency, and, for L ^ S, Mp does not depend on 5; if L ^ S, 
Mp oc (mi/m2) 3> 1 which could therefore generate a very large number of precessions. In Table Ull we show the 
number of precession cycles Mp and the value of the precession frequency fJp at fa and /isco for a number of choices 
of the mass of a binary system. 

Using Eq. H32|l we define now the precession angle a as 



da 



= % , (40) 



From Eqs. ()31|l . and (|35|l it is straightforward to derive the relevant equations for the evolution of the total and 
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orbital angular momentum: 



L + TS 



(1 + 2kT + T2)1/2 ' 
J = L{l + 2T + T'^f'^ , 



= T 



J 
t 

L = JqcosAl 



S(1 + ET)-L(k-T) 
(1 + 2kT + T2)3/2 ' 

Jo X L + e(Jo X L) X L 
z — Jo cos 6 J 



„ sinA^sina 

— — sm Al cos a + (Jo x z) ■. — 

sm Uj sm 0j 

In the previous equation 6j and (pj are the polar coordinates of Jo in the fixed frame and 

Al = arcsm— 



(41a) 

(41b) 
(41c) 

(41d) 

(41e) 

(42) 



is the angle between L and j (at the first order in e, the approximation in which we are working, A^ coincides with 
the angle between L and Jq) and 

S 



T{t) 



From Eqs. H32|) and (|34|l one can write explicitly: 



a,^(2 + |!^) 

Zmi ) r-^ 



(43) 



(44) 



16 /m\3/2 



5 V r / [\ Ami ) 

where 

g^lyj) EE 1 + 2kT + T2 
Eqs. lO, HU, (EHIl and (gHl) can be used to integrate Eq. lB?Hl : 

96 3m2 



1 



S/i'^m-^ \ 4mi 
■'3kS'^ (l - K^) arcsinh 



2(^^i)^ -3KS'(i + K5) 
L + kS 



5(1-^2) 



1/2 



(45) 



(46) 



(47) 



In Eq. H47|l the angle ac is a constant of integration that essentially identifies the position of L and S on the precession 
cone at the reference time tc- Finally, the frequency evolution of Al - the angle between L and J changes with time 
- is given by: 



sin A; 



cos A; 



T 1 - k2 



(1 + 2kT + T2)1/2 

1 + kT 
(1 + 2kT + T2)1/2 



S (1-^2)1/2 

L{f) + nS 

L{f)Q{f) ■ 



(48a) 
(48b) 



Notice that at the 1.5PN order, which corresponds to the the post- Newtonian order at which we model the waveform 
in this paper, the contribution of the spins to the GW phase H23|l depends only on the parameter fi. However, if 
one takes into account the change of orientation of the angular momenta, then the spin-orbit coupling is described 
by three parameters, say k and S, on which /3 depends. In the implementation of the computation of the Fisher 
information matrix, cf Sec.0 we found convenient to use as independent parameters P, k and ac- 
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IV. THE OBSERVED SIGNAL 



We can now derive ready-to-use analytical expressions for the signal measured by LISA at the Michelson detector 
output, Eq. 0, when the source is an in-spiralling binary system of spinning massive black holes. In this case, due 
to the relativistic spin-orbit coupling, the binary angular momenta, L, Si and S2, precess, following the equations 
described in Sec lIII Bl It is worth spelling out the assumption under which the signal is derived: (i) the binary is 
in circular orbit; (ii) the waveform is modelled using the restricted post^'^-Newtonian approximation (amplitude at 
the lowest quadrupole Newtonian contribution, and GW phase at the 1.5PN order); (iii) the masses and/or spins are 
such that the system undergoes " simple precession" . 

In the most general case the signal depends on 17 parameters: 2 mass parameters, 6 parameters related to the BH 
spins - the magnitude and orientation of each spin - the orbital eccentricity, the luminosity distance (which is a direct 
observable in GW astronomy), 2 angles identifying the location of the source in the sky, two angles that describe the 
orientation of the orbital plane, one angle that describes the orientation of the ellipse in the orbital plane, an arbitrary 
reference time, say the time at coalescence, and the signal phase at that time. In the approximation that we consider 
here the signal depends on 12 unknown parameters (cf the next section for more details]. The key new qualitative 
feature that we introduce here and that is not present in any of the previous analysis |2^ |23, HOl is the change 
of orientation of the orbital angular momentum during the observation. Unlike the model considered in [27i l29l | , here 
we consider only the waves emitted at twice the orbital frequency. 

Before providing the explicit expression for h(t), it is worth summarising the key motions and the associated time 
scales that affect the detector output: 

• LISA orbits around the Sun, so that the barycentre of the instrument changes position over the timescale of 
one sidereal year: this motion Doppler-shifts the incoming gravitational waves, and depends only on the source 
position in the sky; 

• The orientation of the LISA arms changes over the same time scale because the detector precesses around the 
normal to the Ecliptic: this introduces a phase and amplitude modulation which is due to the time-dependent 
response of the detector to the two orthogonal polarisations and hx ; this effect depends on both the position 
and orientation of the source in the sky. 

• The binary orbital plane changes orientation in the sky with L and S that precess around the total angular 
momentum J, whose direction is essentially constant. The time scale of precession is longer than the orbital 
period but much shorter than the LISA rotation time-scale and depends on both the source-detector relative 
orientation and the binary physical parameters (masses and spins): the incoming GW's are therefore intrinsically 
modulated in amplitude and phase. 

We consider first the representation of detector output in the time domain /i(')(t); we will then turn to the frequency 
domain representation }S'-\f). 



A. The measured signal in the time domain 



In the amplitude-and- phase representation, the signal h'-''^{t), cf Eq. lO, measured by the t-th Michelson LISA 
interferometer and produced by the polarisation amplitudes (I18all and (jl8b|) reads 



/i(') (t) = A{t) 4') (t) cos U(0 + 6pcj){t) + (t) + (t) 



In the previous expression 




1+ L-N 



- \ 2 



1/2 





l.n) 


F'x\t) 


1 + 


(l.n; 


)1 


> 



are the polarisation amplitude and phase, respectively; 

ipD (t) = 27ri?e / sin On cos($ (t) ~ 4>n) (i?© = 1 AU) 



(49) 

(50) 
(51) 

(52) 
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FIG. 1: The evolution of the polarisation amplitude (left panels) and phase (right panels), computed according to Eqs. H5()|l 
and l|51|l . for the final year of in-spiral of a binary system as a function of time. The bold solid-line refers to black holes with no 
spin, the dotted-line (for the phase) and thin solid-line (for the amplitude) to spinning objects. The relevant source parameters 
are chosen as follows: mi = l(f Mq, m2 = l(f Mq, S/rn^ — 0.95 and S • L = 0.5 (top panels), and mi = m2 — 10^ M©, 
S/m? = 0.3, and S • L = 0.9 (bottom panels). The location and initial orientation of the source have been chosen randomly. 
See the text for further details. 



is the Doppler phase modulation induced by the motion of the LISA detector around the Sun; 5p4>{t) is the so-called 
Thomas precession phase, Eq. H36|l . which vanishes for spin-less binary systems; the amplitude of the gravitational 
wave signal is 

A[t{f)]=2'-^[iTf{t)f'\ (53) 

and the GW phase (/)(t) is given by Ea. (|23|l . where at the 1.5PN order one neglects the term proportional to (ttM/)*/'^. 

The structure of h^'-\t) clearly shows that the signal measured at the detector output is phase and amplitude 
modulated because of the three effects that we have mentioned above: (ii) the orbital motion of the instrument 
around the Sun which Doppler-modulates, through Lpo, the phase of the incoming gravitational wave signal; this 
effect is the same for h'^^^ and ft,'^^^-'; (ii) the change of orientation of the detector during the time of observation 
that through F^^\t) and F^^\t) affects A^'^ and '^^p \ in different ways depending on t; (iii) the change of L due to 

spin-orbit interaction which affects A^''^ and tpp^ through (N • L) and ■f/'^, cf. Eqs. l|5a|) . ((8H)l . ifT^ . ISUJl and (|FT|l : the 
impinging signal is intrinsically modulated. Furthermore 5p4i adds to the phase of the signal. 

In order to complete the analysis of the detector output, we need to provide explicit expressions for (L • N), which 
enters Eq. H5Q(I and (|51(l . and (L • z') and [N • (L x z')] that are needed to compute -0^, Eq. (|12|l . and therefore the 
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antenna beam patters. Using Eqs. I@J and (|41e|l we obtain: 

f ivT KT^ \ , cos 6*^ - (Jo • N) cos 6*7 . 
L-N = (Jo-N)cosAiH ■. — smALCOsa 

sin W 7 



L z' 



sin 6 J 

, (Jo X z) ■ N . . 

H : — sm A L sm a , 

sm Oj 

n -'^ X , l~2(Jo-zOcosgj . , 

= (Jo • z ) cos A/, H smALCOsa 

2 sm Oj 



(Jo X z) • z' . 
: — sm A L sm a , 



(54) 



(55) 



sm 9 J 

<r iCt -/^ , N • (z X z') - N • (Jo X z')cos6ij . ^ 

N-(Lxz') = N • (Jo X z')cosAl H ^ . ^ ^smAiCOsa 

sm Oj 

N • (Jo X z) X z' . ^ . 

H sm A L sm a . 

sm Oj 

sinAi, cosAi, and are a are given by Eqs. H48a|l . (|48b|l and H47|l . respectively, and the 
ition of the source with respect of the detector read: 

(58) 
(59) 



(56) 
(57) 



In the previous expressions sinA^, cosAi, and are a 
remaining quantities describing the orientation of the 

1 \/3 
Jo-z' = -cos 0.7 ^ sin6'jcos($(i) - 0j) , 

Jo • N = cos 0j cos + sinflj sin6'Ar cos(0j — 0jv) , 



N.(. 



Jo X z' 



= - sin 6'Ar sin 6*7 sin(0,7 — 0jv) 



V3 

( 

2 

V3 



\/3 

cos $(t) (cos 6',/ sin sin^jy ~ cosOjy sinOj sin <j>j) 
On sin Oj cos j — cos Oj sin ^ cos ) , 



- — sin$(t) (cos^AT 



N • (Jq X z) = sin sin 0j sin((/)j — 
\/3 

N • (z X z') = — sin0Ar sin [$(t) - 4>m] 
\/3 

J • (z X z') = — sin0j sin [^{t) - 0j] , 



N • (Jo X z) X z 



Lastly, we need to compute the 
Eqs. and P7 |) : 

{Lx N)-L^Vl. 



-— smO J 
2 



V3cos0jv cos [$(t) — 0j] + sin0jv cos(0,7 — 0jv 
scalar and vector products necessary 

cos Al (N • L) - (N • Jo) - e Jo • (L X N) 



to derive the Thomas precession 



(60) 
(61) 
(62) 

(63) 
(64) 

phase Sp(j>{t), 
(65) 



T r-f KT^ N-(JoXz) COS0,7(N • Jo) - COS0Ar . . 

Jo ■ (L X N) = — smAi cosa H — smA^ sma ; 



sm Oj 



sin Oj 



(66) 



Sp4>{t) can be therefore computed using Eqs. (gT)), (|48a|l . (|48b|) . and (|^ . 

An example of the LISA output when black holes rapidly spin is shown in Figure ^ where we plot the amplitude 
and phase modulation, Eqs (|5nil and (|5T|l . during the final year of in-spiral. Two cases are actually presented: a 
binary system with mi = 10^ M©, m2 — 10^ M©, S/m? = 0.95 and S • L = 0.5 and an equal mass binary where the 
parameters correspond to mi = m2 = 10^ M©, S/w? = 0.3, and S • L = 0.9. The evolution of Ap{t) and (/fp(i) for 
the two sets of physical parameters is compared to the spin-less case. A simple eye inspection shows how dramatic is 
the effect induced by spins, and it should therefore not come as a surprise the fact that the errors associated to the 
measurements of the source parameters are strongly affected. 
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B. The measured signal in the frequency domain 

For several applications in data analysis, it is often more useful to work in the frequency domain, and we will derive 
now an approximation to the Fourier representation of the signal h^''\t), Ea. (|49|l . Our convention for the Fourier 
transform g{f) of any real function g(t) is: 

/oo 
e'^^f'g{t)dt. (67) 
-oo 

The Fourier transform H67I) of the signal ft,'^')(i) can be computed in an rather straightforward way using the 
stationary phase approximation, see (52i] . and considering the gravitational waveform A(t) cos (l)(t) in Eq. (|49(l as the 
carrier signal modulated by the motion of the detector and the source's orbital plane l3.l| . Under this assumption, 
the Fourier transform of the detector output H49|) reads: 

hi^)(^f) ^ I AAi^\t{f)] f-y^ ei^ < / < /isco ^ (68) 

jo / > /isco 



where 



*(/) = - 0, - J + J {SttM/)-^/' 



^ 20 /743 11 \ , ^,2/3 

IH \ ri] {■KrafY'^ 

9 V336 47^ 

n^ , r^ ,„/3058673 5429 617 , \, ^x4/3 

4 (47r - B) iTimf) + 10 \ 77 H »? - cr (vrm/)^/^ 

^ ^' ^ ■" V 1016064 1008 ' 144 ' ^' 



(70) 



Notice that in Eq. H70() the term proportional to {-Kmf)'^/'^ correspond to the post^-Newtonian order, and will not be 
considered in this work. 



V. PARAMETER ESTIMATION 



In this section we discuss the errors associated to the parameter measurement when LISA monitors binary systems 
of rapidly rotating massive black holes in circular orbit. We start by briefiv recalling the general concepts and formulae 
regarding parameter estimation ~ we refer the reader to [43, l49l ISOl ISlL l53 | and references therein for more details - 
and then present the results of our analysis applied to LISA data. 



A. Review of signal analysis and parameter estimation 

The signal s{t) registered at the detector output is the superposition of noise n{t) and gravitational waves h(t; A): 

s{t) = h{t;\) + n{t)- (71) 

A represents the vector of the unknown parameters (location, masses, spins, etc.) that characterise the actual wave- 
form and that one wishes to estimate from the data stream. We assume the noise to be stationary and Gaussian, 
characterised by a noise spectral density 5'„(/). In the geometrical approach to signal processing it is useful to consider 
sit) as a vector in the signal vector space, and to introduce the following inner product between two signals v and 
w |53 |: 

According to the definition ifT^ . the optimal signal-to- noise ratio (SNR) at which h can be detected is 



(73) 
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In this paper we discuss the errors associated to the measurement of the unknown parameter vector A that characterises 
the signal h{t;X). In the hmit of large SNR, which is clearly the case for LISA observations of massive black hole 
binary systems, the errors AA follow a Gaussian probability distribution: 



p{AX) = 



det(r) 
2tt 



1/2 



In Ea. (|74|) the matrix Tjk is known as the Fisher information matrix, and reads 



-(0 ^ 

- jk — 



The variance- covariance matrix is simply given by the inverse of the Fisher information matrix: 

jk 



(74) 



(75) 



(76) 



The matrix S contains full information about the parameter errors and their correlations, and is what we need to 
compute (cf next Section) in order to investigate the LISA parameter estimation. In fact the diagonal elements of S 
represent the expected mean squared errors 



((AA^) 



(77) 



and its off-diagonal elements provide information about the correlations among different parameters through the 
correlation coefficients c'^: 



kk 



(-1 < c^^ < +1) , 



(78) 



In the limit of high signal-to- noise ratio, S^-^ provides a tight lower bound to the minimum mean-squared error 
((AA-')^) the so-called Cramer-Rao bound |43,|53- It is important to notice that the errors (|77|l and the correlation 
coefficients (|78ll depend both on the actual value of the signal parameter vector A. For the case of observations with 
two or more detectors with uncorrelated noise, the Fisher information matrix is simply: Tjk — X), ^j^- 

One of the properties that we are interested in is the angular resolution of the instrument, which we define as: 



AilN = 27r |(Acos6'^) (A0^) - (A cos6'jv A0Ar)^| 



1/2 



(79) 



The physical meaning of Afl^ is the following: the probability of the source to lie outside an (appropriately shaped) 
error ellipse enclosing a solid angle Afl is simply ^-'^^/'^^n ^ 



B. LISA observations 



We discuss now how accurately LISA can measure the source parameters in the case of rapidly spinning black holes. 
The goal of this section is to investigate whether spins, which have been neglected in studies carried out so far, can 
significantly affect the estimation of the source parameters. A thorough exploration of this effect requires to probe a 
very large multi-dimensional parameter space, and it is well outside of the scope of this paper. Here we concentrate 
on a fiducial source of two 10^ Mq black holes at redshift z = 1 - a rather typical LISA source - and we compare the 
errors associated to the parameter measurements in the case in which spins are present - for large spins, S/M'^ = 0.95, 
and moderate spins S/M"^ — 0.3; in both cases k = 0.9 - and the case where spins are neglected (the black holes are 
considered, a priory, to be not spinning). 

Before presenting the results wc spell out the assumptions under which we compute the errors |(77J): 

• The orbit is circular; we regard such hypothesis as realistic, as we are dealing with massive systems of black holes 
of comparable mass. As the orbit shrinks due to dynamical friction the eccentricity is likely to decrease [llllT^IT^ 
and radiation reaction completes the circularisation process before G W's enter the observational window of LISA; 
in fact the eccentricity e evolves according to e oc J^is/is [4^^ 
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The black holes are spinning, so that the binary orbit precesses in space. However, in oder to simplify the 
description of the precession motion, still addressing a realistic astrophysical scenario, we assume that either 
nil = (the case considered in this paper) or 5*2 = and the spins and angular momentum are not anti-aligned. 
Under this conditions a binary undergoes the so-called simple precession 33] , and the equations describing the 
evolution of the relevant physical quantities simplify considerably (cf Sec. IIIIBl for more details). 

We restrict the analysis only to the in-spiral phase of the whole coalescence. We approximate the waveform 
h'-''\f;X) at the restricted 1.5PN order. The GW signal detected at the LISA output is therefore described 
by Eqs. H68|l . H69() . and (|70() . Because we retain terms up to 1.5PN order in the GW phase, in the expression 
of ^'(/), Eq. H70|l . we neglect the last term proportional to (ttM/)'*/'^. We shut-off the in-spiral waveform at 
the frequency /isco) given by Eq. (|15|l . In the computation of the errors (|77|l we actually neglect the Thomas 
precession phase Sp(t), Eq. H36|) . This simplification is motived by computational reasons, and does not affect in 
any significant way the final results. In fact, Sp[t{f)] must be computed numerically at each frequency using the 
past history of the binary; 6p[t{f)] needs then to be included in the integrand of the scalar product f 17^ that 
leads to determination of the elements of the Fisher information matrix. This double numerical integration, 
which needs to be carried out to high accuracy in order to keep under control numerical instabilities that 
occur otherwise in the numerical inversion of Fjfe, makes the computational time very long, due to our limited 
computational resources. We have checked for a few (random) choices of the source parameters that including 
Sp[t{f)] does not change in any appreciable way the result. This is to be expected. For the physical parameters 
considered in this paper the Thomas precession phase contribute to « 1 wave cycle, out of a total « 1500. It 
also represents a secular increase in the phase of the GW signal, qualitatively not different from the one given 
by *(/). 

We consider the fiducial sources to be at redshift z = 1 [s^. As the systems are at cosmological distance, 
the values of all the physical parameters entering the GW signal presented in the previous sections must be 
considered as the observed ones; they differ from the values of the parameters as measured in the source rest 
frame by a factor (1 + z); the parameters are " Doppler-shifted" according to: 

/ 

(1 + ^)' 
(l + z)i, 

il + z)M, 

(1 + z) m , 

{l + z)^l. (80) 



/ 



t 

M 

m 



We assume that the instrument observes the whole final year of in-spiral; the GW's sweep therefore the frequency 
band between the arrival frequency and the final cut-off frequency /isco; /a is determined so that after 1 yr, 
as measured by an observer in our solar system, the GW instantaneous frequency reaches /isco: 

^obs ^(/isco) ^(/a) 

= lyr = 3.1556926 X lO'^ sec (81) 

The range /a < / < /isco determines the integration domain in H75|) ; values of /a and /isco for selected choices 
of the source parameters are given in Table |l| Note that we do not impose a low-frequency cut-off to LISA; this 
choice is based on the fact that the real low frequency noise wall of space-based instruments is not very well 
understood at the moment, placed somewhere in the range 10~^ Hz - 10""* Hz. 

The total noise that affects the observations is given by the superposition of instrumental sources and astro- 
physical foregrounds of unresolved radiation due to (mainly) galactic white dwarf binary systems 53, 5 J| , the 
so-called "confusion noise". The total noise spectral dei isity Sn{f) is therefore the sum of these two components, 
and we use the analytical approximations given in |26l l55| . 

Out of the 17 parameters on which the most general waveform depends, the signal h^^\f\ A) that we consider 
here depends on 12 independent parameters. In our analysis we adopt the following choice of independent 
parameters: InA^ and In/i (mass parameters), cosOn, <Pn, cosOj, (pj, and ac (geometry of the binary with 
respect to the detector), n and f3 (spin parameters), InD^ (distance of the source) and finally tc and (f>c- 

1/2 

We compute the expected mean square errors ((AA^)) and the angular resolution Ail^ for the case of 
observations carried out with one detector and with combined observations of both detectors / and //. The 
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FIG. 2: The probability distribution of AQn (left) and ADl/Dl (right), for observations of the final year of in-spiral of a 
binary system with mi — m2 = lO'' Mq at z = 1. The histograms show the result of a Monte-Carlo simulation, where 1000 
sources have been randomly located and oriented in the sky. The bottom panels refer to measurements carried out with only 
one detector, whereas the top panels describe the results obtained by combining the two independent data streams. The plots 
compare the errors for different values of the BH spins: S/m^ — 0.9 (solid line), 0.3 (dotted line), and (dotted-dashed line). 



analysis is done in the the frequency domain: we first compute analytically the derivatives 9/i(''/5A^, where 
j — 1, .., 12, then compute numerically Tjk and E^*^, Eqs. H75|l and (|76|l : the integration and matrix inversion 
are performed using numerical routines of the NAG library. 

LISA parameter estimation strongly depends (the results vary by orders of magnitude) on the actual value of the 
signal parameters. In particular, one of the key set of parameters that affect the errors is the location and orientation 
of a source with respect to the detector. This represents already a large parameter space that one needs to explore 
in order to obtain meaningful results. We perform this exploration by means of Monte-Carlo simulations, which 
therefore affect the way in which the results are presented, given in terms of probability distributions. For each set 
of physical parameters, we select randomly the 5 geometrical parameters {On, (j)^, 9j, (pj and ac) from a uniform 
distribution in cosfl^v, 0jv, cos9j, and ac- The Monte-Carlo simulation is done on 1000 different sets of angles. 
As far as the other 7 parameters are concerned we chose them as follows. We consider a fiducial source at redshift 
z = 1, with nil — 1712 = 10^ Mq (which sets the 3 parameters Dl, A4 and /it), with tc — (pc — 0. We fix the tilt angle 
parameter k — 0.9, and we explore three different values of the size of the spin S: (i) S/M"^ = 0.95, (ii) S/M"^ = 0.3 
and (iii) 5' = 0. Notice that given mi, m2, S and n we can easily derive (3. It is worth noticing that in this paper we 
do not explore the effect of the tilt angle cos~^ k, which is also likely to affect the errors. Such analysis is currently 
in progress 0. We would also like to stress that for the case 5 = the waveform that we actuall y c onsider is the 
one corresponding to the non-precessing restricted 2PN approximation (this is the waveform used in [33| , very similar 
to the one used in |26j |. which is the non-precessing restricted 1.5PN approximation). This is equivalent to assuming 
that we know a priori that spins are zero and the binary system does no precess. Note that this is different from 
estimating the source parameters using the waveform (|68f) where precession is included, and assuming S <^1. In the 
case 5 = we need to estimate only 11 parameters, and not 12. The reason of this choice is to compare existing 
results in the literature with our new ones that take into account spin-orbit modulations. 

Figs m and 21 summarise the key results. Figs |3 13 show the probability distribution of the parameter errors, for 
Ar^AT, H^DlI Dl, M. and A/x//x. Fig ^ shows the cumulative probability distributions: if we define ^ = AA^, 

the plots show P(^ < ^o) = /q^" v{^)^^ as a function of ^o- Each plot compares the non-spinning case with the cases 
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FIG. 3: The probability distribution of AAi/M (left) and A/j./ /i (right) for observations carried out with only one interferometer 
(bottom) and the two combined LISA outputs (top). The parameters are the same as in Fig|2| 



in which the black holes are spinning with two different value of the spins, S/M'^ — 0.95 and 0.3 . 

The key result, which is absolutely clear from Figs |21 El a-nd ^ is that the errors in the spinning-case are smaller 
than in the non-spinning one, even if the number of parameters that one needs to estimate is greater (12 instead of 
11). One can intuitively understand this behaviour by looking at Figure ^ which shows the amplitude and phase 
evolution of the signal recorded at the LISA output: precessing binaries show a much greater richness of features than 
non-spinning sources. These are the features that help in measuring the parameters. 

However, despite all the parameters are estimated with smaller errors with respect to the non-spinning case, there 
is striking difference between the "position parameters" - location in the sky and luminosity distance - and mass 
parameters. In fact, the errors AJItv and /^Di^/D^ are reduced by a factor 2-to-lO depending on the actual value 
of the source parameter. We know that LISA reconstruct the position of the source in the sky by exploiting the 
modulation in the amplitude and phase of the signal, and precession introduces one additional modulation effect 
which can improve the determination of N and de-correlates Dj^ from M, which improves the measurement of the 
luminosity distance. Still, the change of position and orientation of LISA is the main effect that allows the source 
to be located. The improvement in Ailjv and ADl/Dl is however significant: the fraction of sources that can be 
located within, say, one square degree is 10% for non-spinning binaries, 35% for binaries with 0.3 and 50% with 0.95 
(for K = 0.9). The systems whose distance is known to better than 1% is about 60% for non-spinning binaries, and 
essentially the totality for spinning binaries. 

The situation is radically different for the masses: in the non-spinning case the instrument reconstructs the value 
of the two mass parameters by staying in phase with the GW phase, whose time evolution is mainly controlled by M 
(the leading Newtonian term) and, to a less extent, by /i, through the post-Newtonian corrections, cf Eqs. (|23|l . H70|) 
and Table m M and play no role in the amplitude and phase modulation, Ap{t) and (pp{t), Eqs l|50|) an d 1511). As 
a consequence, Ai is measured much more accurately, by a factor ^ 100 or more, than fi, see Fig ^ and |2a,|33- If 
spins are present, the masses start to play a role also in Ap{t) and (pp{t), because they control the rate at which L 
and S precess, and therefore are responsible for the intrinsic amplitude and phase modulation of the signal. This has 
two effects on the parameter estimation: it de-correlates Ai from fi (at the Newtonian order the two mass parameters 
are even degenerate) and provides new features in the measured signal that improve the parameter estimation. The 
net effect is that the errors in M and fj, are now of the same order, with a drastic improvement with respect to the 
non-spinning case, of by a factor ~ 50 and ~ 10'^ for AM/Ai and Afi/fj,, respectively; AM/Ai is still 3 times 
smaller than A/i/^, due to the role played by Ai in the GW phase (j>{t). It also clear that the additional signature 
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FIG. 4; Same as Figure 2, but now the cumulative probability distributions as a function of the parameters AQn and ADl/Dl 
are shown. See text for more details 



produced by the masses on the amplitude and phase modulation helps in removing the correlation between M and 
Dl, therefore reducing ADl/Dl, as we have already mentioned. 

To summarise, the parameters of two 10^ M© black holes spiralling toward the final merger at 2; = 1, and rapidly 
rotating, say S/M'^ = 0.95 and k = 0.9, can be measured very accurately: A4 and /i can be measured to a few 
parts in 10^, the luminosity distance to better than 1% (in some cases almost 0.1%), and with an error box in the 
sky 10~^ srad ;S AJl^v ~ 10~^srad. Notice that the errors scale with the distance roughly as Ailjv ~ and 
AAi/A4 ~ An/fj, ~ ADl/Dl ^ ^/D^. This would be strictly true for white noise, which is not the case for LISA. 
In particular going to higher redshift, one can expect a degradation of the measurements more severe than the one 
predicted by this simple scaling, because GW's are red-shifted to lower frequencies where the noise is higher and some 
initial portion of the signal falls out of the observational window. 

For multi-band, electro-magnetic and gravitational, observations of the same event the key parameter is the instru- 
ment's angular resolution. In fact there is some hope that LISA will have enough angular resolution to locate the 
galaxy or galaxy cluster where a coalescence of massive black holes takes place; this would enable other telescopes to 
be pointed at the same area of the sky and to observe the aftermath of such catastrophic event. However, previous 
studies l2ik ,30 J have concluded that this might actually be impossible for sources at z « 1 (or beyond), 

because of the poor angular resolution of GW observations. In short, the results of our analysis suggest that in some 
exceptional case LISA might be able to identify the host galaxy cluster of a MBH binary of lO^M© at z = 1, but in 
general the LISA error box will contain several hundreds of galaxy clusters and 10^ galaxies, all of them potential 
hosts of a LISA detected source. Such conclusion can be easily derived by comparing the LISA angular resolution 
with the one of other telescopes and the angular size of galaxies at high redshift. An optical telescope, such as the 
Keck telescope, has a field of view of about 1 square degree. Chandra has a square field of view of 16 x 16 arcmin with 
a spatial resolution of ~ 1 arcsec, and XMM has a circular field of view with radius 15 arcmin and a spatial resolution 
of a few arcsec. INTEGRAL has an even wider circular field of view - the radius is 10 degree - with a poorer spatial 
resolution (12 arcmin). There is therefore no doubt that, for a typical LISA source at z = 1, GW observations will 
provide information accurate enough to cover the interesting portion of the sky with a single observation by one of (or 
all) the former instruments. At z = 1, the typical size of a galaxy is ~ 3 arcsec and the typical size of a galaxy cluster 
is ~ 2 arcmin. In a very few lucky occasions LISA will be able to detect a source with AQ,n w 3 x 10^^ deg^, but more 
typically the angular resolution is AQn ~ 0.3 deg^, and could be much worse for a few detections (Ail^v ~ 5deg^). It 
is therefore clear that for the "high spatial resolution" GW detections, only « 3 galaxy clusters will fall in the LISA 
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error box, but more typically the number of clusters will be « 300. If there is still a small (but not negligible) chance 
of identifying the galaxy cluster hosting a GW source, the chances of pin-pointing the host galaxy seem to be very 
bleak (the former numbers need to be multiplied by a factor lO'^ for galaxies). 

VI. CONCLUSION 

We have considered LISA observations of massive black hole binary systems in the final stage of in-spiral for 
rapidly rotating black holes. We have restricted our analysis to comparable mass objects in circular orbit, modelling 
the radiation at the restricted post -'^■^ -Newtonian order assuming simple precession. We have derived read-to-use 
analytical expressions for the signal registered at the detector output, both in the time and frequency domain, and 
have determined the mean-square-errors associated to the parameter measurements for equal mass lO^M© binary 
systems at z = 1 for selected spin parameters and a wide range of source locations and orientations. Our analysis 
clearly shows that the presence of spins reduces (by orders of magnitude, for some parameters) the errors with which 
the source parameters are measured. LISA is therefore a more powerful telescope than previously thought if spins 
play a significant role. The main shortcoming of our analysis is the limited region of the parameter space that we 
have been able to explore, which is entirely due to our limited computational resources. Such analysis is currently 
in progress [3^ . Of particular interest is the exploration of the parameter space in the case where the masses are 
not comparable, i.e. the mass ratio is 0.1 or smaller. Our present understanding of the relevant astrophysical 
scenarios suggests, in fact, that the formation of unequal mass binaries should be regarded as the rule rather than the 
exception. When the mass ratio increases, the number of precession cycles detected at the LISA output will increase 
as well, cf. Tabic Hll and Eq. H39|) . It is therefore conceivable that precession will be even more effective in breaking 
the degeneracy among the parameters, in particular the two masses. However the amount of intrinsic rotation S shall 
also play a crucial role, as the number of precession cycles is proportional to 5'. Preliminary results suggest that for 
TO2/TO1 = 0.1 and rapidly rotating systems {S/M^ « 0.9) the estimation of the mass parameters and the angular 
resolution improve by a factor of a few, with respect to the equal mass case, provided that the signal-to-noise ratio is 
not strongly affected. In general, in fact, unequal mass systems will be detected at a lower SNR, so that precession 
compensates for the loss of SNR. The latter eventually dominates the effect of rotation when m2 /mi decreases further 
and compromises the accuracy with which parameters can be measured. 
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